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ABSTRACT 

We compute mass outflow rate R m from relativistic matter accreting quasi-spherically 
onto Schwarzschild black holes. Taking the pair-plasma pressure mediated shock sur- 
face as the effective boundary layer (of the black hole) from where bulk of the outflow 
is assumed to be generated, computation of this rate is done using combinations of 
exact transonic inflow and outflow solutions. We find that R m depends on the initial 
parameters of the flow, the polytropic index of matter, the degree of compression of 
matter near the shock surface and on the location of the shock surface itself. We thus 
not only study the variation of the mass outflow rate as a function of various physical 
parameters governing the problem but also provide a sufficiently plausible estimation 
of this rate. 
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1 INTRODUCTION 

Computation of mass outflow rates from the advective 
accretion disks around black holes and neutron stars very re- 
cently has been done (Chakrabarti, 1998; Das, 1998, 1998a; 
Das & Chakrabarti, 1998 [DC98 hereafter]), by self consis- 
tently combining the exact transonic accretion and wind 
solutions. Rigoriously justifying the fact that most of the 
outflowing matter comes out from the CENtrifugal Pressure 
Supported BOundary Layer (CENBOL) it has been shown 
there that the Rankine-Hugoniot shock location or the loca- 
tion of the maximum polytropic pressure acts as the CEN- 
BOL. However, for some black hole models of active galactic 
nuclei, inflow may not have accretion disk (Rees, 1977 and 
references therein) . Accretion is then quasi-spherical having 
almost zero or negligible angular momentum (Bondi [1952] 
type accretion) and the shock is not of Rankine-Hugoniot 
type. On the otherhand, absence of angular momentum rules 
out the possibilty of formation of the polytropic pressure 
maxima. So for quasi-spherical accretion, absence of intrin- 
sic angular momentum of the accreting material does not al- 
low CENBOL formation. It has been shown that (Meszaros 



& Ostriker, 1983; Kazanas & Ellison, 1986 [KE86 hereafter]) 
for quasi -spherical accretion onto black holes, steady state 
situation may be developed where a standing collisionless 
shock may form due to the plasma instabilities and for non- 
linearity introduced by small density perturbation. This is 
because, after crossing the sonic point the infalling matter 
(in plasma form) becomes highly supersonic. Any small per- 
turbation and slowing down of the infall velocity will cre- 
ate a piston and produce a shock. A spherically symmetric 
shock produced in such a way will accelerate a fraction of 
the inflowing plasma to relativistic energies. The shock ac- 
celereted relativistic particles suffer essentially no Compton 
loss and are assumed to lose energy only through p — p col- 
lision. These relativistic hadrons are not readily captured 
by the black hole (Protheros & Kazanas, 1983) rather con- 
sidereble high energy density of these relativistic protons 
would be maintained to support a standing, collisionless, 
spherical shock around the black hole (see KE86 and ref- 
erences therein). Thus a self-supported standing shock may 
be produced even for accretion with zero angular momen- 
tum. In this work, we take this pair-plasma pressure medi- 
ated shock surface as the alternative of the CENBOL which 
can be treated as the effective physical hard surface which, 
in principle mimics the ordinary stellar surface regarding 
the mass outflow. The condition necessery for the develop- 
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ment and maintainence of such a self supported spherical 
shock is satisfied for the high Mach number solutions (Elli- 
son & Eichler, 1984). Keeping this in the back of our mind, 
for our present work, we concentrate only on low energy 
accretion to obtain high shock Mach number. Considering 
low energy (£ £ 0.001) accretion, we assume that particles 
accreting toward black hole are shock accelerated via first 
order Fermi acceleration (Axford, et. all, 1977) producing 
relativistic protons. Those relativistic protons usually scat- 
ter several times before being captured by the black hole. 
These energised particles, in turn, provide sufficiently out- 
ward pressure to support a standing, collisionless shock. A 
fraction of the energy flux of infalling matter is assumed to 
be converted into radiation at the shock standoff distance 
through hadronic (p — p) collision and mesonic (^,71°) de- 
cay. Luminosity produced by this fraction is used to obtain 
the shock location for the present work. Our approach to for- 
mulate the expression for shock location is somewhat similar 
to that of KE86. 

At the shock surface, density of the post-shock material 
shoots up and velocity falls down, infalling matter starts 
piling up on the shock surface. The post shock relativistic 
hadronic pressure then gives a kick to the piled up mat- 
ter the result of which is the ejection of outflow from the 
shock surface. For this type of inflow, accretion is known to 
proceed smoothly after a shock transition, since successful 
subsonic solutions have been constructed (Flammang, 1982) 
for accretion onto black holes embedded within normal stars 
with the boundary condition u = c; where u is the infall ve- 
locity of matter and c is the velocity of light in vaccum. 
The fraction of energy converted, the shock compression ra- 
tio Rcomp (in the notation of DC98), along with the ratio 
of post shock relativistic hadronic pressure to infalling ram 
pressure at a given shock location are obtained from the 
steady state shock solution of Ellision and Eichler (Ellison 
& Eichler, 1984; 1985). The shock location as a function of 
the specific energy £ of the infalling matter and accretion 
rate is then self consistently obtained using the above men- 
tioned quantities. We then calculate the amount of mass 
outflow rate Rm from the shock surface using combination 
of exact transonic inflow outflow solutions and study the de- 
pendence of Rfn on various physical entities governing the 
inflow-outflow system. We thus quantitatively compute the 
mass outflow rate only from the inflow parameters. In this 
way we analytically connect the accretion and wind type 
topologies self-consistently. As there was no such attempt 
available in the litereture which computes the mass loss rate 
from zero angular momentum quasi-spherical Bondi type ac- 
cretion, our work, for the first time we believe, could shed 
light on the nature of the outflow from the models of AGNs 
with no accretion disks. 

We consider polytropic inflow. The outflow is also assumed 
to be polytropic except the fact that ~f out fi ow is assumed to 
be less than the yinfiow reason of which will be discussed 
in the next section. As a fraction of the energy of infalling 
material is converted into radiation, energy flux of the wind 
is somewhat less than that of the accretion but is kept con- 
stant througout the outflow. 
The plan of the paper is as follows : 

In the next Section, we describe our model and present the 
governing equations for the inflow and the outflow along 
with the solution procedure. In §3, we present the result of 



our computation. Finally, in §4, we draw our conclusions. 
A preliminary report of this work has been presented else- 
where (Das, 1998b). 



MODEL DESCRIPTION, GOVERNING 
EQUATIONS AND THE SOLUTION 
PROCEDURE 



2.1 Inflow model 

We assume that a Schwarzschild type black hole of mass 
IOM0 quasi-spherically accretes low energy (£ 0.001) 
fluid obeying polytropic equation of state. The density of 
the fluid is p(r), r being the radial distance measured in the 
unit of Schwarzschild radius. We also assume that the ac- 
cretion rate (in the unit of Eddington rate) with which the 
fluid is being accreted, is not a function of r. For simplicity 
of calculation, we choose geometric unit (unit of length = 
2l ^ M & unit of velocity is c. M is the mass of the black hole, 
c is the velocity of light and G is the universal gravitational 
constant. All other relevant physical quantities can be ex- 
pressed like wise) to measure all the relevant quantities. We 
ignore the self gravity of the flow and the calculation is being 
done using Paczyn'ski-Wiita (Paczyn'ski & Wiita, 1980) po- 
tential which mimics surronding of the Schwarzschild black 
hole. The equations (in dimentionless geometric unit) gov- 
erning the inflow are : 

a) Conservation of specific energy is given by: 



£ = 



u(r) 



+ na(r) — 



(1) 



2(r-l) 

u(r) and a(r) are the radial and polytropic sound veloci- 



ties respectively. a(r) = 



p(r) 



p{r) being the polytropic 



pressure. For a polytropic inflow, p(r) — Kpir) 1 where K is 
a measure of the entropy which remains constant throught 
the flow except at the shock location, n is the polytropic 
constant of the inflow n = (7 — 1) _1 , 7 being the polytropic 
index. 

b) Mass conservation equation is given by, 

2 



M in = <d in p(r)u(r)r 



(2) 



As already mentioned, we assume that a steady, collision- 
less shock forms at a distance r sh (measured in the unit of 
Schwarzschild radius) due to the instabilities in the plasma 
flow. We also assume that for our model, the effective thick- 
ness of the shock A sh is small enough compared to the shock 
standoff distance, i, e, 

Ash « r sh 

Accreting particles with infall velocity u(r) are then assumed 
to be shock accelerated via first order Fermi acceleration. 
Due to this process, relativistic protons will be produced. 
These relativistic protons, suffer essentially no Compton 
loss and hence are not readily swallowed by the black hole. 
Rather they usually scatter several times before being cap- 
tured by the black hole thus provide sufficient outward pres- 
sure necessary to support the shock. These protons, in turn, 
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produce pions through post-shock inelastic nuclear collisions 

P + P- 



± i o 
p + p + n + tv 



Pions generated by this process, decay into relativistic elec- 
trons, neutrinos & antinutrinos and produces high energy 7 
rays. 



27 



+ v e 



e + + u e 



These electrons produce the observed non-thermal radiation 
by Synchrotron and inverse Compton scattering. The overall 
efficiency of this mechanism depends largly on the shock 
location. It has been shown (Eichler, 1979) that almost half 
of the energy flux that goes into relativistic particles is lost 
due to neutrinos. 

From eqs. (2), the density of infalling matter p(r) comes out 
to be 



p(r) 



Qi„u(r)r 2 



(3) 



At the shock, density of matter will shoot up and inflow 
velocity fall down abruptly. If (p_,u_) and (p + ,u+) are the 
pre and post-shock densities and velocities respectively, then 



— Re 



(4) 



Where R CO mp is the shock compression ratio (in the nota- 
tion of DC98). For high shock much number solution (which 
is compatible with our low energy accretion model), the ex- 
pression for R com p can be well approximated as 



Rcomp = l.UM 3h ' 



(5) 



Where M s h is the shock Mach number and eqs. (5) holds 
good for M sh £ 4.0 (Ellison & Eichler, 1985). However, the 
scattering mean free path of the relativistic hardrons pro- 
duced by this process are assumed to be small enough so 
that they could encounter the full shock compression ratio 
while crossing the shock. 

The hadronic interaction charecteristic time scale t pp may 
be expressed as 

Tpp = — (6) 

p+CJppC 

Where a pp is the collision cross section for relativistic pro- 
tons. Using eqs. (4), t pp can also be expressed as 



1 



Rcomp P — 0~pp C 



(6a) 



* As only a fraction of the accreting matter is shock enerzyzed, 
the value of p+ used to calculate t vv is, in reality, less than that 
of the actual p+ , giving a higher value of t pp . The accurate value 
of Tpp can be calculated using the cosmic ray energy spectrum 
and coupling the relativistic and non-relativistic part of the ac- 
creting plasma which will be presented elsewhere. Nevertheless, 
our rough estimation assures that even with the accurate value 
of Tpp (which is higher than that of used here), the conditions of 
shock formation are highly satisfied. 



If we now assume that a fraction £f of the infalling energy is 
converted into radiation through the hadronic collision (p - 
p) and mesonic {-k^ 1 ,tt°) decay, this £f will allow convergent 
steady state solutions and in the case of quasi-spherical in- 
fall, allows the development and maintenance of a standing, 
collisionless shock at a fixed distance r s h measured in the 
unit of Scwarzschild radius. 

The luminosity obtained from this energy £f at the shock 
location r s h is assumed to be L\. For a spherical shock sur- 
face and taking care of neutrino losses, the expression for Ci 
can be written as (KE86) 

(7a) 



£1 = 2nr 2 a hPshm p u 3 s h£F 



where m p is the mass of the proton and subscript sh indi- 
cates that the respective quatities are measured at the shock 
location. 

If we assume that the pressure of the relativistic particles 
Vrci (uniform inside the shock), the average energy density is 
then 3Vrci , so alternatively the luminosity can be expressed 
in terms of the volume integral of the the emissivity e due 
to the hadronic (p — p) collition where 

3Vrel 

e = 

Tpp 

Thus the alternative expresion for the luminosity obtained 
as function of 7-Vei would be, 



(76) 



Defining 5 as the ratio of downstream relativistic particle 
pressure to incoming ram pressure at the shock, we obtain 



Prel = 8p s hm p u 



(8) 



Equating (7a) and (7b) and putting the values of p, Vrel , 
and t pp from eqs. (3), (6a) and (8) respectively, we get the 
expression for the shock location as a function of various 
inflow parameters as 



_ 2>a vv cMi n ( 6 

U s h 2 Qin \£f 



(9) 



The ratio (^r-J as a function of shock Mach number M s h 
for a high shock Mach number solution (low energy inflow) 
is obtained from the empirical solution deduced by Ellison 
and Eichler (Ellison & Eichler, 1984), after suitable modifi- 
cation required for our model. 



2.2 Outflow model 

In ordinary stellar mass loss computations (Tarafder, 1988 
and references therein) the outflow is assumed to be isother- 
mal till the sonic point. This assumption is probably justi- 
fied, since copious photons from the stellar atmosphere de- 
posit momenta on the slowly outgoing and expanding out- 
flow and possibly make the flow close to the isothermal. This 
need not be the case for outflow from black hole candidates. 
Our effective boundary layer, being pretty close to the black 
hole, are very hot and most of the photons emitted may 
be swallowed by the black hole itself instead of coming out 
of the region and depositing momentum onto the outflow. 
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Thus, the outflow could be cooler than the isothermal flow 
in our case. We choose polytropic outflow with a different 
polytropic index j < 7 due to momentum deposition. As a 
fraction of infalling energy density (£f) is converted into ra- 
diation, specific energy of the outflow is somewhat less than 
that of the inflow. Nevertheless, the outflow specific energy 
is also kept constant throughout the flow. 
The following two conservation laws are valid for the outflow 



lution: 



c' v ( r ) 2 , ' 1 \i 
£ — — h n a(r) — 



1 



2(r - 1) 



M out = Qoutp{r)v(r)r 



(10) 



(11) 



Where £ ' is the specific energy of the outflow and £ ' < £ and 
n' = (70 — l) -1 is the polytropic constant of the outflow. 
&out is the solid angle subtended by the outflow and v(r) is 
the velocity of the outflow. 

For simplicity of calculation, we assume that the outflow is 
also quasi-spherical and & ou t~Oin- Defining Rrh as the mass 
outflow rate, we obtain 



Rri 



M 



M ln 



(12) 



It is obvious from the above discussion that Rrh should have 
some complicated functional dependences on the following 
parameters 

Rrh = *(£ , M ln ,r sh , Rcom P r/, 7°) (12a) 



2.3 Procedure to solve the inflow and outflow 
equations simultaniously 

Before we proceed into further detail, a general understand- 
ing of the transonic inflow outflow system in present case 
is essential to understand the basic scheme of the solution 
procedure. Let us consider the transonic accretion first. In- 
falling matter becomes supersonic after crossing a saddle 
type sonic point, location of which is determined by the in- 
flow parameters such as the specific energy £ , Mi n (in the 
unit of Eddington rate) and 7 of the inflow. This supersonic 
flow then encounters a shock (if present) location of which 
(r s h) is determined from eqs. 9. At the shock, part of the 
incoming matter, having higher entropy density, is likely to 
return back as wind through a sonic point other than the one 
through which it just entered. Thus a combination of tran- 
sonic topologies, one for the inflow and other for the outflow 
(passing through a different sonic point and following the 
topology completely different than that of the "self-wind" 
of the accretion), is required to obtain a full solution. So 
it turns out that finding a complete set of self-consistent 
inflow outflow solutions ultimately boils down to locating 
the sonic point of the outflow and the mass flux through it. 
Thus a supply of parameters £,M;„,7 and 70 make a self 
consistent computation of Rrh possible. Here y is supplied 
as free parameter because the self-consistent computation of 
7o directly using £ , M;„ and 7 has not been attempted in 
this work, instead we put a constrain that 7 < 7 always 
and for any value of 7. In reality, y is directly related to 
the heating and cooling procedure of the outflow. 
The following procedure is adopted to obtain a complete so- 



From eqs. (1) & (2), we get the expression for the deriva- 



tive 



du(r) 
dr 



2(r-l)^ 



u(r) y4- 

V ' u{r) 



(13) 



At the sonic point, the numerator and denominator sepa- 
rately vanish and give rise to the so called sonic point con- 
dition; 



u c 



2(r c - 1) 



(14) 



Where the subscripts c represents the quantities at the sonic 
point. The derivative at the sonic point (^r) c is computed 
using the L' Hospital rule. The expression for is ob- 

tained by solving the following polynomial, 



' 2n + 1 \ / du \ 2 3u c / du \ 

/\drj c nr c \drj c | (r c — l) 3 



2a e 2 / 



71+ 1 

n 
(15) 

Using fourth order Runge Kutta method, tt(r) and a(r) 
are computed along the inflow from the inflow sonic point 
till the position where the shock forms. The shock location 
is calculated by simultaniously solving the eqs. (1), (2) and 

(9) . With the known value of £ ' and j a , it is easy to com- 
pute the location of the sonic point of the outflow from eqs. 

(10) and (11). At the outflow sonic point, the outflow ve- 
locity v c and polytropic sound velocity a c is computed in 
the same manner performed for the inflow. Using eqs. (10) 
and (11), (^) and (^) is computed as was done for the 
inflow. Runge-Kutta method is then employed to integrate 
from the outflow sonic point towards the black hole to find 
out the outflow velocity v and density p at the shock loca- 
tion. The outflow rate is then computed using eqs. (12). 



3 RESULTS 

Fig. 1 shows a typical solution which combines the accretion 
and the outflow. The input parameters are £ — 0.001, M n 
= 1.0 Eddington rate (£d stands for the Eddington rate in 
the figure) and 7 = | corresponding to relativistic inflow. 
The solid curve with an arrow represents the pre-shock re- 
gion of the inflow and the solid vertical line with double 
arrow at X pps (the subscript pps stands for pair plasma me- 
diated shock) represents the shock transition. Location of 
shock is obtained using the eqs. 4 for a particular set if in- 
flow parameters mentioned above. Three dotted curves show 
the three different outflow branches corresponding to differ- 
ent polytropic index of the outflow as 7 = 1.3(left most 
curve), 1.275 (middle curve) and 1.25(rightmost curve). It 
is evident from the figure that the outflow moves along the 
solution curves completely different from that of the "wind 
solution" (solid line marked with an outword directed arrow) 
of the inflow which passes through the sonic point P a . The 
mass loss ratio Rrh for these cases are 0.0023, 0.00065 and 
0.00014 respectively. In Fig. 2., we have plotted the variation 
of Rrh with incoming specific energy £ for a set of values of 
M in (measured in the unit of Eddington rate shown as £ d in 
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the Fig.) shown in the Fig. It is observed that R m monoton- 
icaly increases with energy. This is because as £ increases 
keeping the Eddington rate of the inflow fixed, the shock 
Mach number M sh decreases result of which is the decrea- 
ment of shock location r s h and post shock density (via eqs. 
5.) but the increament of the post shock fluid velocity (v sh ) 
with which the matter leaves the shock surface. The out- 
fow rate R m , which is the product of these three quantities, 
in general increases monotonically with £ due to the com- 
bined tug of war of these three quantities. Morever, closer 
the shock forms to the black hole, the greater will be the 
amount of gravitational potential available to be put onto 
the relativistic hadrons to provide more outward pressure 
at the shock boundary which gives a stronger "kick" to the 
accreting matter, the result of which is the increament in 
Rm- All these points are manifested in Fig. 3 where we have 
shown the variation of R m as a function of compression ratio 
Rcomp (solid curve), the shock location r s h (dotted curve) 
and the injection velocity of the outflow v 3 h (dashed curve) . 
The fig. is drawn for a fixed 7 = | and -y = 1.3 Rcomp and 
v rsh are scaled as R comp — > {Rcomp — 5.890) x 10 3 and v rsh 
4 x W~ 6 v rs h- The unequal gaps between the curves with 
different £a implies that when inflow energy £ is kept fixed, 
Rm nonlinearly increases with Eddington rate. This is be- 
cause, as £ is kept fixed while Mj„ is varied, the amount of 
infalling energy converted to produce high energy protons is 
also fixed, so higher is the value of M in (in the unit of Ed- 
dington rate), the larger is the distance of the shock surface 
(measured from the black hole) and the outflowing matter 
feels low inward gravitational pull, the result of which is the 
non- linear corelation of R m with M in . 

To have better insight of the behavior of the outflow, in Fig. 
4a. and 4b, we plot R m as a function of the polytropic index 
of the inflow 7 (Fig. 4a) and that of the outflow 70 (Fig. 4b) 
for fixed £ = 0.001 and M in — 1.0£d- The range of 7 shown 
here are the range for which shock forms for the specified £ 
and Mi n . In Fig. 4b, it is observed that shock forms for the 
values of 7 lower than that shown in the figure, but mass 
loss rate is then so small that we did not plot for 7 < 1.23. 
The general conclusion is that R m correlates with 7 . This 
is because as 70 increases, shock location and post shock 
density of matter does not change (as 7 does not have any 
role in shock formation or in determining the Rcomp) but 
the sonic point of the outflow is pushed inward, hence the 
velocity with which outflow leaves the shock surface goes up 
resulting the increament in R m - However, R m anticorelates 
with 7 which is observed from Fig. 4a. 



4 CONCLUDING REMARKS 

In the present paper, we have computed the mass outflow 
rate from the relativistis low energy matter quasi-spherically 
accreting onto Schwarzschild type black holes. The free pa- 
rameters chosen for the inflow are the specific energy £ , 
mass inflow rate Mi„ (In the unit of Eddington rate) and 
the polytropic index 7 of the inflow. Only one extra parame- 
ter which was supplied for the outflow is its polytropic index 
7o- We have computed the mass outflow rate R m taking the 
help of inflow parameters only (except 7 ) thus could an- 
alytically connect the accretion and wind type topologies 



self-consistently. In our computation, we could investigate 
the dependence of R m as a function of different physical 
quantities governing the inflow. 

The basic conclusions of this paper are the followings, 

(i) It is possible that outflows for quasi-spherical Bondi 
type accretion onto a Schwarzschild black hole are coming 
from the pair plasma pressure mediated shock surface. 

(ii) The outflow rate monotonically increases with the 
specific energy of the inflow and nonlinearly increases with 
the Eddington rate of the infalling matter. 

(iii) R m , in general, corelates with 7 but anticorelates 
with 7. 

(iv) Generaly speaking, as our model deals with high 
shock Mach number (low energy accretion) solutions, out- 
flows in our work always generate from the supersonic 
branch of the inflow, i,e, shock is always located inside the 
sonic point. 

(v) Unlike the mass outflow from the accretion disks 
around black holes (Das, 1998; 1998a; DC98) here we found 
that the value of R m is distinguishably small. This is because 
matter is ejected out due to the pressure of the relativis- 
tic plasma pairs which is less enough in comparison to the 
pressure generated due to the presence of significant angu- 
lar momentum. However, in the present work we have dealt 
only high Mach number solution which means matter is ac- 
creting with very low energy (cold inflow, as it is described 
in litereture). This is another possible reason to obtain a low 
mass loss rate. If, instead of high Mach number solution, we 
would use low Mach number solution, e,g, high energy ac- 
cretion, the mass outflow would be considerably higher (this 
is obvious because it has already been established in present 
work that R m increases with £ . In our next work, we will 
present this type of model by calculating R CO mp and (^7) 
for low Mach number solution. 

In the literature, we did not find any numerical simulation 
work which deal with the type of outflow we discuss in this 
paper (outflow from zero angular momentum inflow). On the 
other hand, observationally it is very difficult to calculate the 
outflow rate exactly from a real system, as it depends on too 
many uncertainities, such as filling factors and projection ef- 
fects. So, at present, we don't have any offhand results in 
this field with which our result could be compared. 

There are a number of possible improvements which 
could be made on this preliminary work. For instance, the 
effect of radiation pressure on both the inflow and the out- 
flow could be taken into account. A preliminary investiga- 
tion shows that the effect of radiation force, when included 
into the basic conservation equations (eqs. (1), (2) and (10), 
(11)) decreases the value of R m - This is probably because 
introduction of any radiation term which is proportional to 
Jj, wekens gravity and pushes the shock location outwards, 
combined effect of which is the decreament of R m - Another 
possible improvement is to include the magnetic field to give 
the outflow an appropiate geometry. In our model, we as- 
sumed the outflow to be quasi-spherical as like the inflow. 
The introduction of magnetic field probably could collimate 
the outflow thus would help to investigate the structure of 
jet related to the model of AGN without accretion disk in 
greater detail. Lastly, we did not self consistently compute 
7o as a function of inflow parameters and y was supplied as 
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a free parameter. In our future work, we will be presenting 
the self consistent calculation of y to reduce the number of 
free parameters of the problem. 

So far, we made the computation around Schwarzschild 
black hole. Our work could be extended to study mass 
outflow in kerr space time using pseudo-Kerr potential 
(Chakrabarti & Khanna, 1992) and incorporating the frame 
dragging effect. This is under preparation and will be pre- 
sented elsewhere. 
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Figure Captions 



Fig. 1: Solution topology for three different y (1.3, 1.275, 
1.25) for £ = 0.001, M in = 1.0, £ d = 1.0, 7 = f . P s indicates 
the sonic point of the inflow where X pps stands for the shock 
location. See text for details. 

Fig. 2: Variation of R m with inflow specific energy £ for dif- 
ferent Eddington rate (marked as £d in the figure). Other 
parameters are 7 = | and 7 =1.3. R m monotonically in- 
creases with £ , see text for detail. 

Fig. 3: Variation of Rm with the compresion ratio R CO m P 
(solid curve), shock location X ps (dotted curve) and outflow 
velocity at shock v s h (dashed curve). Other parameters are 
£=0.001, M in = 1.0(In the unit of Eddington rate), 7 = § 
and 7 =1.3. 
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Fig. 4a: Variation of R m with polytropic index of inflow 7. 
Other parameters are £=0.001, M in = l.Q (in the unit of Ed- 
dington rate) and 70 =1.3. 

Fig. 4b: Variation of R m with polytropic index of outflow 
7o. Other parameters are £=0.001, Mi n =1.0 (in the unit of 
Eddington rate) and 7 = |. 
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